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We have searched for proton decay via p — > p + K° using data from a 91.7 kiloton-year exposure of Super- 
Kamiokande-I, a 49.2 kiloton-year exposure of Super- Kamiokande-II, and a 31.9 kiloton-year exposure of 
Super-Kamiokande-III. The number of candidate events in the data was consistent with the atmospheric neu- 
trino background expectation and no evidence for proton decay in this mode was found. We set a partial lifetime 
lower limit of 1.6x 10 33 years at the 90% confidence level. 
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I. INTRODUCTION 



The experimental observation of proton decay would pro- 
vide strong evidence for Grand Unified Theories (GUTs) [1] 
which unify the strong, electromagnetic, and weak forces. Su- 
persymmetry (SUSY) GUTs predict the unification of cou- 
pling constants consistent with the non-observation of p — > 
via gauge boson exchange. In some SUSY GUT 
models [3] p — > vK + is then favored by dimension five op- 
erators. However, in other models [4] p — >• p + K° may dom- 
inate. So far there has been no experimental observation of 
either of these decays and in this paper we focus on the later, 
whose most stringent partial lifetime lower limit comes from 
Super-Kamiokande (SK), z/B^+go > 1.3 x 10 33 years (90% 
C.L.) |Ht]. The partial lifetime of p —> p + K° is predicted to be 
just above this limit 0]. 

The K° is a composite state of the K$ (50%) and K% (50%). 
The analysis presented in [5] searched only for proton decay 
into p + K® using data from the first run period of SK (91.7 
kton-years). In this paper we present a combined search for 
proton decay into both K$ and with an updated data set 
(172.8 kton-years). 

A p — > p + K° signal would appear in SK as one muon- 
like (/j-like) ring with monochromatic momentum of about 
327 MeV/c and secondary ring(s) from the decay of the K°. 
The total momenta of all such rings should be consistent with 
the decay of a proton, as well as an invariant mass close to 
the proton mass. The dominant decay modes of the K% are 
JI+7T (69.2%), and 7t°7t° (30.7%). The K% decays predomi- 
nantly into i&e^Me (40.6%), nV^V (27.0%), 37t° (19.5%), 
and TC+TTTI (12.5%). Since the K® has a relatively long decay 
length (cpx ^8.4 m, where x is the mean life), proton decay 
events into this channel are characterized by two vertices sep- 
arated in time and space: a primary vertex from which the p + 
emerges and a secondary vertex for the decay products of the 
K®. We developed new event reconstruction for events with 
two vertices. 

The outline of this paper is as follows. A review of the SK 
detector is presented in Section II. In section III the p — > p. + K° 
Monte Carlo (MC) simulation as well as the atmospheric neu- 
trino MC simulation used to calculate the signal efficiency and 
expected background rate are presented. Section IV and V 
describe the data reduction and the event reconstruction, re- 
spectively. The results of the data analysis are summarized in 
Section VI and a comparison with the previous result |5J is 
described in subsection VI-C. Finally, Section VII concludes 
this study. 
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II. SUPER-KAMIOKANDE DETECTOR 

The Super-Kamiokande detector [6] is located in Gifu 
Prefecture, Japan, at a depth of 2,700 m water equivalent 
within the Kamioka mine at Mt. Ikenoyama. The detector 
is a 50 kton right cylinder that is optically separated into 
two concentric regions, an inner (ID) and an outer detector 
(OD). Inward-facing 20 inch diameter Hamamatsu photomul- 
tiplier tubes (PMTs) 0] are mounted uniformly on the walls 
of the ID. During the SK-I (SK-III) run period from April 
1996 (October 2006) to July 2001 (September 2008) approx- 
imately 11,100 PMTs lined the ID walls. There were about 
5,200 PMTs in the ID during the SK-II run period from Octo- 
ber 2002 to October 2005. The photocathode coverage during 
the run was about 40% (20%) in SK-I and -III (SK-II). A black 
sheet spans the region between PMTs to absorb light and op- 
tically separate ID and OD. Starting with SK-II, each ID PMT 
is encased in an acrylic cover to reduce the effect of implosion 
on neighboring PMTs. The transparency of this cover is more 
than 96% for >350 nm photons with normal incidence. A 
22.5 kton fiducial volume used in the data analysis presented 
below is defined by the cylindrical volume taken 2 m from the 
ID walls. 

The OD is a 2 m thick cylindrical shell enclosing the ID 
that is instrumented with 1,885 outward-facing 8 inch diam- 
eter Hamamatsu PMTs. To improve the light collection effi- 
ciency, a 60 cm x 60 cm wavelength shifting plate is attached 
to each PMT and the OD walls are covered with reflective 
material called "Tyvek sheet". The main purpose of the OD 
is to tag incoming cosmic ray muons and muons induced by 
atmospheric neutrinos that escape the ID. The OD region also 
serves as a passive shield against radioactivity from outside 
the detector wall. The SK-I detector and its calibration are 
described in detail elsewhere |@]. 



III. SIMULATION 

The p -» p + K° MC and the SK atmospheric neutrino MC 
are used to estimate signal efficiencies and the expected back- 
ground rates. Using these MCs the event selections described 
in Sec. VI-A are optimized to maximize the analysis' sensi- 
tivity prior to examining the data. 



A. Proton Decay 

The p — > p + K° MC used here is the same as that used in 
05J with the exception that p — > p + K [ l events are addition- 
ally simulated in this analysis. In this analysis the decays of 
both free proton and protons bound on the oxygen nucleus are 
considered. For the decay of free protons the directions of 
the /j + and K° are exactly back-to-back, and each has a mo- 
mentum of 326.5 MeV/c. In the case of proton decay within 
the oxygen nucleus, the effects of Fermi momentum, corre- 
lation with other nucleons, the nuclear binding energy, kaon- 
nucleon interactions, and K® — > K9 regeneration are all taken 
into account. The Fermi momentum distribution measured by 
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electron- 12 C scattering H] is used in the simulation. The nu- 
clear binding energy effect is taken into account by modifica- 
tion of the proton's mass. The modified proton mass, m', is 
calculated by m' p — m p — Eb where m p is the proton rest mass 
and Eb is the nuclear binding energy. The value of Eb for 
each simulated event was randomly selected from a Gaussian 
probability density function with (p, o) = (39.0, 10.2) MeV 
for the S-state and (ji, a) = (15.5, 3.82) MeV for the P-state. 
The correlation of the decaying proton's wave function with 
other nucleons reduces the invariant mass of the decay prod- 
ucts since some of the system's momentum is carried by these 
nucleons. This effect is present in about 10% of proton de- 
cays (9J and produces a broad invariant mass spectrum down 
to about 600 MeV/c 2 . The location of protons within the lfi O 
nucleus is calculated according to the Woods-Saxon nuclear 
density model II 1(111 . The kaon nucleon interactions within the 
16 nucleus that are considered include elastic scattering and 
inelastic scattering via charge exchange. The K + N scatter- 
ing amplitudes were fitted by a partial wave analysis using 
many data samples by Hyslop et al. Hill . The K + n — !• K°p 
charge exchange cross section was measured by Glasser et 
al. B12I1 . The detail of the K + N interaction is described in 
01 • From isospin symmetry, the K°N reactions have essen- 
tially the same magnitude as the K + N reactions. If an inter- 
action occurs, the effect of Pauli blocking is introduced by re- 
quiring the nucleon momentum after interaction be larger than 
the Fermi surface momentum, given by pf{t) — (j7t 2 p(r))3, 
where p (r) is the nuclear density distribution and r is the dis- 
tance from the center of the nucleus. The — > K2 regen- 
eration effect is implemented as a decay mode of the K® in 
this analysis. The regeneration probability in water is based 
on the results of a kaon scattering experiment using a carbon 
target [13]. Regeneration within the ie O nucleus is also con- 
sidered and is assumed to be proportional to its density relative 
to water. The fraction with regeneration is about 0. 1 % of the 
total generated p — > P + K® MC events. 

We simulate propagation of particles and Cherenkov 
light production in water by custom software based on 
GEANT 01411 . The propagation of charged pions in water is 
simulated using custom code based on [15] for pion momenta 
less than 500 MeV/c and using CALOR 1 16] otherwise. 

B. Atmospheric Neutrinos 

The interactions of atmospheric neutrinos within SK pose 
the dominant source of background to searches for nucleon 
decay. For the p — > p + K° mode, the decay products of the 
K® and K® often include pions which can be created in atmo- 
spheric neutrino interactions. Notably, charged current mul- 
tiple pion production and deep inelastic scattering processes, 
again with multiple pions in the final state, are backgrounds 
to the search presented here. Accordingly, the accurate sim- 
ulation of atmospheric neutrino events in SK is essential to 
properly estimate the proton decay background. In this paper, 
the atmospheric neutrino MC is the same as used in other SK 
analyses j p , UJn . The primary flux is taken from the Honda 
model [18] and atmospheric neutrino interactions are gener- 



ated within the detector using the NEUT simulator lfl9i l20ll . 
Kaon regeneration is not considered in the atmospheric neu- 
trino MC because the fraction of K° production in these events 
is negligibly small. For each of the SK run periods presented 
here, independent about 200 year MC samples were used to 
compute the background rates for the present search. Two- 
flavor neutrino oscillation effects with sin 2 2 923 = 1.0 and 
Am 2 = 2.5 x 10~ 3 eV 2 are included in the background predic- 
tion. A more detailed explanation of the atmospheric neutrino 
MC is presented in lUfnll. 

IV. DATA SET AND REDUCTION 

The current data set corresponds to a 91.7 ktonyears 
(1489.2 days) exposure from SK-I, a 49.2 kton years 
(798.6 days) exposure during SK-II, and a 31.9 kton years 
(518.1 days) exposure from SK-III. The data acquisition trig- 
ger threshold used in this data analysis during these periods 
corresponds to the energy deposited by 6 MeV electron dur- 
ing SK-I and SK-III and 8 MeV electron during SK-II. Most 
events above this threshold are cosmic ray muons or low en- 
ergy backgrounds from radioactivity in detector materials. To 
remove these backgrounds the data are passed through sev- 
eral reduction stages before detailed reconstruction processes 
(described in Sec. V) are applied. To summarize, the starting 
point of the data analysis is the fully contained (FC) fiducial 
volume (FV) data sample, which contains events passing the 
following cuts: 

• the number of hit PMTs in the highest charge cluster of 
PMTs in the OD is less than 10 

• the total visible energy in the ID is greater than 30 MeV 

• the reconstructed vertex's distance to the ID wall is 
greater than 2 m 

After these cuts the FCFV event rate is about eight events per 
day. Non-atmospheric neutrino backgrounds are estimated to 
represent less than 1% of this sample. Details of the data re- 
duction can be found in lfl7ll . 

V. EVENT RECONSTRUCTION 

Event reconstruction programs are applied to the FCFV 
sample described in the previous section. For the K® searches, 
the standard event reconstruction IT7i l2lll used in the SK 
atmospheric neutrino oscillation analyses and proton decay 
searches is used. However, relative to the previous version of 
this analysis 05J], several improvements to the reconstruction 
algorithms have been made. The search for K® events makes 
use of new event reconstruction methods, which are explained 
in detail below. 

The event reconstruction for the searches is: 

• vertex reconstruction 

The event vertex is initially determined by using the 
time of flight (TOF)-subtracted timing distribution of 
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hit ID PMTs assuming all light originated from a sin- 
gle point (point fit). An initial direction for the event 
is then estimated as the sum of vectors drawn from this 
vertex to each hit PMT and weighted by its observed 
charge. Defining this direction as the particle's trajec- 
tory, the distribution of observed charge as a function of 
the opening angle to this axis is used to determine the 
edge of the brightest Cherenkov ring in the event. With 
the ring edge in hand, the vertex position and particle 
direction are refined using the timing residual distribu- 
tion, but assuming that PMTs within the ring are illumi- 
nated by photons generated along the particle track and 
other PMTs receive light generated at the event vertex. 

• ring counting 

Based on the reconstructed vertex, additional ring can- 
didates in the event are searched for using a technique 
based on the Hough transform [ 22] . Possible ring candi- 
dates returned by the Hough mapping are tested using a 
likelihood method to determine if the existence of a sec- 
ond ring better describes the observed PMT hit pattern 
than a single ring alone. If a second ring is found, the 
procedure is iterated up to a total of five rings. At each 
step the likelihood is based on a comparison of the ob- 
served charge in the event against the expectation from 
the number of proposed rings. With respect to the previ- 
ous analysis there have been improvements to the prob- 
ability density function (PDF) for the expected charge 
that are used in the present analysis. More details are 
presented below. 

• PID (e//j separation) 

After the number of rings has been determined each 
ring's particle type (PID) is estimated taking into ac- 
count the effects of light scattering in the water and re- 
flection off PMT surfaces. Two likelihood functions are 
used to separate electromagnetic shower type (e-like) 
rings from muon type (/j-like) rings. The first is a pat- 
tern matching algorithm that exploits topological dif- 
ferences in the PMT hit pattern that arise from the scat- 
tering and pair production processes an electron expe- 
riences, but that the primarily unperturbed muon does 
not. A second method makes use of the difference in 
the observed Cherenkov angle between the two parti- 
cles as a result of their largely different masses. For the 
final estimation of a ring's PID, the two methods are 
combined. 

• decay electron search 

Electrons produced from the decay of muons in the de- 
tector can occur in the same trigger timing gate with the 
primary event or delayed. The search for decay elec- 
trons extends to both cases. 

• momentum reconstruction 

The momentum of each ring is determined using the 
sum of its assigned charges after correcting for the at- 
tenuation of Cherenkov light in water and the incident 



angle dependence of the PMT acceptance. Depending 
upon the result of the PID estimation, the momentum is 
then assigned based on the amount of corrected charge 
expected for the fitted particle type. A pionic momen- 
tum (ft-like) that relies on the expected Cherenkov an- 
gle as well as the total corrected charge is also used in 
this analysis. The time variation of the light attenuation 
length and the PMT gain, as measured with cosmic rays 
traversing the detector, is incorporated into the momen- 
tum reconstruction. 

Improvements made to the single photoelectron response 
functions used to predict the expected charge in each PMT 
have led to improvements in the ring counting, PID, and mo- 
mentum estimations for this analysis. 

The uncertainty in the detector's energy calibration (the en- 
ergy scale) is one of the most important for the proton de- 
cay searches. It is estimated by systematic studies of the 
differences between the absolute momenta of data and MC 
in several control samples. At SK these control samples are 
chosen to probe a variety of energies where neutrino interac- 
tions and nucleon decay searches are relevant, among which 
are cosmic ray muons that stop in the detector, their asso- 
ciated decay electrons (Fig. [TJ, and 7t°s produced in neutral 
current atmospheric neutrino interactions (Fig. [2j». The uncer- 
tainty of the absolute energy scale is estimated to be less than 
0.74%, 1.60%, and 2.08% for SK-I, SK-II, and SK-III, re- 
spectively. The time variation of the reconstructed momenta 
of these muon and electron samples during each of the SK 
run periods (Fig. is also taken into account. Overall, the 
uncertainty on the energy scale is 1.1%, 1.7%, and 2.7% in 
SK-I, SK-II, and SK-III, respectively. Further, the directional 
dependence of decay electron momenta is used to estimate 
the asymmetry of the energy scale. During SK-I, SK-II, and 
SK-III, the asymmetry was estimated to be 0.6%, 0.6%, and 
1.3% respectively. These errors on the SK energy scale are 
incorporated into the proton decay analysis as discussed in 
Sec. VI-B-2. 



A. K[ event reconstruction 

The K® event reconstruction needs two vertices, one for the 
fi + candidate from the proton decay and another for particles 
produced in the delayed decay. The difference between the 
two vertices is used to distinguish the K® decay signal from 
the atmospheric neutrino background. However, the standard 
reconstruction described cannot by itself reconstruct two ver- 
tices. For instance, the standard ring counting algorithm uses a 
cut on the PMT hit timing after correcting for the time of flight 
taken from the initial vertex to help isolate genuine Cherenkov 
rings. Since the light from the particles emerging from the 
is both delayed in time and spatially separated from this ini- 
tial vertex, the inclusion of this PMT timing cut degrades the 
performance of the ring counting algorithm on these events. 
Therefore, we developed a specialized reconstruction. In the 
first step of the K® reconstruction, the standard reconstruction 
algorithm is applied without making any cut on the PMT hit 
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FIG. 1: The momentum distribution of decay electrons of the data 
(dot) and MC (solid line) for SK-I, SK-II, and SK-III from top to 
bottom, respectively. MC events are normalized by the number of 
observed data events. The mean values of the data distributions agree 
with those of the MC within 0.6%, 1.6%, and 0.8% for SK-I, SK-II, 
and SK-III, respectively. 



FIG. 2: Invariant mass distribution of neutrino-induced 31 events 
of the observed data (dot) and the atmospheric neutrino MC events 
(solid line) for SK-I, SK-II, and SK-III from top to bottom, respec- 
tively. MC events are normalized by the livetime of the observed 
data. The peak positions of the data distributions agree with those 
from the MC within 0.7%, 1.3%, and 0.3% for SK-I, SK-II, and SK- 
III, respectively. 



timing during the ring counting stage. It outputs an initial ver- 
tex, the number of candidate rings, and a direction for each 
ring in much the same way as the standard reconstruction. 

In the second step of the reconstruction l23ll . an event's hit 
PMTs are divided into two groups for each ring found by 
the first step. The first, called the inner domain, is a col- 
lection of hit PMTs that fall within a cone that is 10° larger 
than the ring's reconstructed Cherenkov angle. The second 
group, called the outer domain, are the remaining collection 
of hit PMTs that are not part of the inner domain. The charge 
and time information of PMTs in the outer domain are zeroed 
and part of the standard reconstruction is reapplied using only 
information from PMTs in the inner domain. At this stage 
the vertex reconstruction, PID, and momentum reconstruction 
are repeated, and the result is stored. A precise vertex fit, 
which adjusts the initial candidate vertex along the ring's re- 
constructed momentum using the expected charge distribution 
for a particle with the ring's new PID, is applied. The vertex 
resolution is about 0.3 m [17]. Finally, an algorithm designed 
to distinguish protons from muons [24], is applied to provide 
an additional PID handle. Once the p. + candidate from the 
original proton's decay is identified, this precision vertex will 
be identified as that of the initial decay. The secondary vertex 
from the decay, is estimated by zeroing the charge and time 
information of the inner domain PMTs and applying the point 
fit algorithm to the outer domain. Since 71 ± from the K% de- 



cay may be just above Cherenkov threshold, these pions may 
not appear in the detector as clear Cherenkov rings. Using the 
outer domain to perform the secondary vertex fit in this way 
incorporates light from even these particles to provide a better 
estimate of the point where the K® decays. These reconstruc- 
tions in the second step are repeated for each ring found by 
the first step. 

For each event, the ring with /j-like momentum closest to 
326.5 MeV/c is chosen as the /j + candidate. A vertex separa- 
tion parameter is defined as the distance between the p. + (pri- 
mary) candidate vertex and the K® (secondary) vertex along 
the /u + candidate's momentum vector. This parameter is used 
to define a cut to help reduce atmospheric neutrino back- 
grounds and is described below. Figure [4] shows the typical 
correlation between the true vertex separation, extracted us- 
ing MC information, and the reconstructed vertex separation 
for p — > p. + K® decay MC events. A clear correlation between 
the two can be seen. A typical event display from one MC 
event showing both the primary and secondary fitted vertices 
appears in Fig. [5] In this case the reconstruction has correctly 
identified the /j + (thick cyan ring) and the primary and sec- 
ondary vertex are represented by the blue triangle and green 
circle, respectively. 

The difference between the reconstructed and true vertex 
separation is shown in Fig. |6]for proton decay events in each 
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FIG. 3: The time variation of the reconstructed momentum/range of 
stopping muons and the reconstructed momentum of decay electrons 
as a function of elapsed days with respect to the mean value for SK- 
I, SK-II, and SK-III from top to bottom, respectively. The largest 
variation (RMS) among these two calibration sources during SK-I, 
SK-II, and SK-III was 0.88%, 0.55%, and 1.79%, respectively. The 
time variation was relatively larger in SK-III due to its poorer water 
quality. 
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FIG. 4: Correlation between true and reconstructed vertex separa- 
tions for the p — > p + K® decay MC events in SK-I after the event 
selections (D1-D5) described in Section VI- A. 




FIG. 5: Typical event display of p — > p. + — ► n 
in SK-I. The upper (lower) visible ring are from true p + (e). The 
cyan thick circle shows the p + candidate ring. The blue triangle and 
green circle show the reconstructed vertex position horizontally and 
vertically projected on the detector wall for the p + and remaining 
particle candidates, respectively. For this event, reconstructed and 
true vertex separations are 6.58 m and 6.33 m, respectively. 



of the SK run periods. A slight bias in this distribution, repre- 
sented by the shift of the peak away from zero, is used to es- 
timate the uncertainty coming from the vertex separation cut 
discussed in Sec. VI-B-2. This shift is typically (9(10) cm for 
both proton decay and atmospheric neutrino MC events. The 
vertex separation resolution for these two MC sets are about 
0.8 m and 1.1m, respectively. Due to the cleaner back-to- 
back event topology of the K% decay events, the separation 
resolution for the proton decay events is slightly better than 
that for the atmospheric neutrino events. Most of this resolu- 
tion comes from the application of the point fit algorithm in 
determining the secondary vertex. 



VI. DATA ANALYSIS 



A. Event Selection 

The search for p — > pi + K° decay is performed by searching 
for decays into /li + K® and /J + K^ separately. The FCFV events 
(Sec. IV) are first passed through the K% selection and failing 
events at the K® selection are then additionally subjected to 
the K® search. The search criterion are not changed with 
respect to the previous analysis @], but are reviewed here for 
completeness. 

To maximize the analysis' sensitivity to the mode, three 
distinct selections designed to extract events corresponding to 
each of the K% decay's final states are used. It should be noted 
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FIG. 6: Reconstructed - true vertex separation for the p — > p + K^ 
decay MC events in SK-I, SK-II, and SK-III (from top to bottom), 
respectively, after the event selections (D1-D5) described in Sec- 
tion VI-A. 



that by construction there is no kinematic overlap among the 
methods. For K® — > n°n° decays (30.7% branching fraction), 
the proton decay selection criteria are as follows: 

• (Al) 3-5 rings, corresponding to the initial p + and the 
gammas emerging from the Jt° decays 

• (A2) one /j-like ring, to be identified as the p + candi- 
date, and the other rings in the event are e-like 

• (A3) one Michel electron, from the decay of the p + 

• (A4) 150 MeV/c < p M < 400 MeV/c 

• (A5) 400 MeV/c 2 < m K o < 600 MeV/c 2 

• (A6) p p < 300 MeV/c 

• (A7) 750 MeV/c 2 < m p < 1000MeV/c 2 

Here p M refers to the p + candidate momentum, is the re- 
constructed K° invariant mass using only the e-like rings, and 
Pp ( m p) 1S the total momentum (invariant mass) using all rings 
in the event. Figure [7] shows the m p and p p distributions af- 
ter applying these criteria (A1-A5). Under the hypothesis of 
proton decay, initial proton momenta of up to 300 MeV/c are 
considered in the search because the Fermi momentum of pro- 
tons in 16 nuclei may reach about 200 MeV/c. This, in addi- 
tion to the possibility of missing gammas from the Jt decay, 
is similarly the reason for the small lower bound on the m p 
cut. 



After this event selection the remaining background con- 
sists primarily of events from v e or Vu charged current (CC) 
multiple pion production interactions, and neutral current 
(NC) deep inelastic scattering processes, from the atmo- 
spheric neutrino MC. These backgrounds often have a single 
/j-like ring back-to-back with several e-like rings and are con- 
sistent with the expected event signature of the signal. 

To search for proton decays in which the K® decays into 
7T, + 7T.~ (69.2% branching fraction) two different methods are 
used. The first method (Method 1) is designed to search for 
events where one of the outgoing pions does not produce a 
clear Cherenkov ring, for instance when its momentum is be- 
low the Cherenkov production threshold, or when it captures 
on a lfi O nucleus. The search criteria are: 

• (Bl) two rings 

• (B2) both rings are /j-like, one from the p + and another 
from one of the charged pions 

• (B3) two Michel electrons, from the decays of the pri- 
mary p + and that from the charged pion 

• (B4) 250 MeV/c < p p < 400 MeV/c 

• (B5) p p < 300 MeV/c 

In this method the most energetic /j-like ring is taken to be 
the p + emerging from the initial proton decay and is assigned 
a /j-like momentum. On the other hand, the second energetic 
ring is assumed to be the charged pion and is assigned a ft-like 
momentum. Figure [8] shows the p p and p p distributions after 
applying these criteria (B1-B3). The second method (Method 
2) searches for proton decays with K® — > n + n~ in which all of 
the emerging charged particles produce Cherenkov radiation 
in the detector. In this case the selection cuts become: 

• (CI) three rings 

• (C2) one or two Michel electrons 

• (C3) 450 MeV/c 2 < m K o < 550MeV/c 2 

• (C4) p p < 300 MeV/c 

• (C5) 750 MeV/c 2 < m p < 1000 MeV/c 2 

Figure [9] shows the m p and p p distributions that result from 
applying the criteria (C1-C3). The primary background to this 
search arises from CC interactions with single- or multiple- 
pions produced in the final state. 

To search for proton decays with a decay the following 
criteria are applied: 

• (Dl) 500(100) pe < total charge < 8000(3000) pe in 

SK-i,m(SK-n), 

• (D2) the number of rings without the PMT hit-timing 
cut > 1 

• (D3) at least one /j-like ring that is neither e-like nor 
proton-like 

• (D4) the number of Michel electrons > 
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FIG. 7: Total invariant mass and momentum for events that satisfy criteria (A1-A5) in SK-I, SK-II, and SK-III from top to bottom. The box 
shows criteria (A6, A7). From left to right, data (1489 days, 799 days, 518 days), atmospheric neutrino MC (200 years, 172 years, 194 years), 
and p — > /j + K° MC are shown. 



• (D5) 260 MeV/c < Pi> < 410 MeV/c 

• (D6) vertex separation between /j + candidate and K% 
decay is greater than 2.3 m 

The /j + -candidate is taken to be the /j-like ring with momen- 
tum closest to 326.5 MeV/c. At the third cut (D3), the follow- 
ing criteria are required to increase purity of the /j-like rings. 



In addition to requiring a /j-like classification by the standard 
PID, the result of the proton ID algorithm must also be /j-like. 
Further, the Cherenkov opening angle is required to be greater 
than 30°, The cut value at (D6) has been chosen to maximize 
this analysis' sensitivity to the proton partial lifetime based 
on studies of the proton decay and atmospheric neutrino MC. 
Figure[l0]shows the vertex separation distribution after apply- 
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FIG. 8: jj + momentum and total momentum for events that satisfy criteria (B1-B3) in SK-I, SK-II, and SK-III from top to bottom. The box 
shows criteria (B4, B5). From left to right, data (1489 days, 799 days, 518 days), atmospheric neutrino MC (200 years, 172 years, 194 years), 
and p -» fj+K° MC are shown. 



ing criteria (D1-D5). The net decay length of the p — > n + K^ 
MC is about 2 m because roughly 70% of the undergo 
hadronic interactions in the water before decaying. There is 
good agreement between the data and the atmospheric neu- 
trino MC. 

The dominant background sources after the event selection 
are primarily CC quasi-elastic and single-pion production 



interactions. In these backgrounds the recoiling proton from 
the neutrino interaction is above Cherenkov threshold and is 
selected as the ^-candidate. For this reason the primary 
vertex can be pulled along the direction of the /j + candidate 
momentum, increasing the separation between the two recon- 
structed vertices. Figure [TT] shows a typical event display of 
these backgrounds. In this particular event the reconstructed 
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FIG. 9: Total invariant mass and total momentum for events that satisfy criteria (C1-C3) in SK-I, SK-II, and SK-III from top to bottom. 
The box shows criteria (C4, C5). From left to right, data (1489 days, 799 days, 518 days), atmospheric neutrino MC (200 years, 172 years, 
194 years), and p -> /j+K° MC are shown. 



vertex separation is 3.9 m. 

B. Result 

The signal detection efficiency, the expected atmospheric 
neutrino background rate, and the number of candidate events 



for each of the search methods are summarized in Table U 

The detection efficiencies in SK-II are lower than in SK-I 
and SK-III and their differences are at most 20% as shown 
in Table I. The estimated background rates of SK-I, SK-II, 
and SK-III are on the order of a few events and the dominant 
source of uncertainty is the MC statistical error as described 
in Sec. VI-B-2. For instance, the statistical error on the back- 
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TABLE I: Summary of the p — > [i + K° search. Systematic errors are shown in the signal efficiencies and background rates (see Table Ull for 
more detail). The numbers in the parentheses of the background columns show the expected background rates in Megaton-years . 




vertex separation (m) 

FIG. 10: Vertex separation distributions for events that satisfy criteria 
(D1-D5) in SK-I, SK-II, and SK-III from top to bottom. Positive 
values of the vertex separation indicate that the point fit vertex is 
separated from the p. + vertex in the direction opposite of the /j + . Data 
(black dot), atmospheric neutrino MC (blue thinner histogram), and 
p — > p + K° MC (red thicker histogram) are shown. The MC samples 
are normalized to data by entries. The cut value corresponding to 
(D6) is shown as broken line. 




FIG. 11: Typical event display of the remaining atmospheric neu- 
trino background MC event in the A*£ search in SK-I. This is from 
charged current quasi-elastic interaction. The right (left) visible rings 
are from true proton (muon). The cyan ring shows the p + candidate 
reconstructed ring. The blue triangle and green circle show the re- 
constructed vertex position horizontally and vertically projected on 
the detector wall for the p + and remaining particle candidates, re- 
spectively. 



ground in the K% — > 71+71 Method 2 search is 20-40%. 

The number of data and the expected background and signal 
efficiencies at each step of the event selection are compared 
for all four searches. For example, Figures [12] shows the K { [ 
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search. After each event selection, the number of remaining 
events in the data is consistent with the background expecta- 
tion in each of the SK run periods. The Poisson probability 
to observe the number of candidates after all event selections 
for each search was calculated. For example, the only proba- 
bility smaller than 5% is seen in the SK-II K® — > 7C°7C° search 
(2%). However, the number of candidates is smaller than the 
expected backgrounds in SK-I and SK-III for this search. The 
sums of the expected backgrounds and the number of candi- 
dates from all SK run periods for the K% and K% searches are 
13.2 and 13, respectively. They are consistent and therefore 
we conclude that there is no evidence for proton decay found 
in this analysis. 

All the candidate events in TableUJwere inspected by hand. 
By checking event display with the reconstructed informa- 
tion, all the events show Sub-GeV (visible energy is less than 
1330 MeV) back-to-back multi visible rings and no obvious 
mis-reconstruction. For example, Figure[l3]shows typical real 
data candidate event in the K® search. This event looks consis- 
tent with the background where the visible proton is iden- 
tified as the /j + candidate ring (see Figure fTTb. 



the background, respectively. Here Eo,; and £>o,i correspond 
to the z'-th searches efficiency and background contamination, 
respectively. These values are summarized in TableQ] The pa- 
rameters a e , and correspond to the total systematic error 
on the signal and background composition for the z'-th search. 
Systematic error estimations are discussed in the next section 
and Table [TT] summarizes the errors for each run period and 
search mode. Finally, A-o,, represents the detector exposure, 
taken to be product of the number of target protons in the de- 
tector and the livetime of each SK run. The corresponding 
uncertainty on the exposure, Ox r is assumed to be 1%. 

We calculate the lower limit of the nucleon decay rate, 
riimit, using a 90% confidence level (CL) as: 



CL = 



/r r =o'nf = if(r|«Q^r 

Sr = oUliP(r\ni)dT 



where N (= 12) is the number of searches. The lower lifetime 
limit: 



x/B 



1 



K° 



Turrit 



1. Lifetime Limit 

In the absence of significant excess of observed data above 
the background expectation we calculate a lower limit on the 
proton partial lifetime using a Bayesian method ll25ll . The cal- 
culation used here follows the one from the previous analysis 

a. 

The decay of the proton is assumed to follow a Poisson 
probability, P(T\nj), and is expressed as: 



p{r\m) = 



^(TXiti+bi) 



{TXfii + bif- 



P{r)P(Xi)P(ei)P(bi)dXideidbi 



where n, is the number of candidate events in the z'-th proton 
decay search (z = 1, 2, 3, and 4 for K% ->■ 7C°7C°, K% — > 7C + 7C~ 
Method 1, K% ->• 7C+7C~ Method 2, and searches, respec- 
tively, in SK-I, z = 5-8 in SK-II, and z = 9-12 in SK-III in Ta- 
ble H|, r is the total decay rate, A,, is the detector exposure, e, 
is the detection efficiency including the meson branching ra- 
tio, and bi is the number of expected background events. The 
probability density function, P(T), describing the prior expec- 
tation of decay rate is taken to be P(T) = 1 for T > and zero 
otherwise. 

Introduction of systematic error effects is done through 
prior probabilities for the detector exposure P(Xi), efficiency 
P(Ei), and background expectation P(bj). They are all as- 
sumed to be truncated Gaussian distributions of the form, 



P( Xi ) • 



(Xi < 0) 

(xi = \i,£i,bi) 



where A-o,/ (Ox ), Eo.i (c*e,) and Z?o,i (Oi,) are the estimates (sys- 
tematic errors) of the exposure, the detection efficiency, and 



is set to be 1.6 x 10 33 years at 90% CL. The result of the limit 
calculation for each search method and data taking period is 
shown in Table U 



2. Systematic Error 



Table [TT] summarizes the systematic errors on the K® 



JC°7E° 



71+71 Method 1 



7C+7IT Method 2. 



and K° 



searches, respectively. 

As for the physics simulation, following uncertainties are 
taken into account. In the signal efficiency, differences in the 
Fermi momentum modeling (8[ [l9ll and uncertainties in the 



correlated decay probability 



the K° L 



K$ regeneration 



probability, as well as the charged pion-nucleon cross section 
in water 15J are considered. The dominant systematic error 
sources are correlated decay and the pion-nucleon cross sec- 
tion. For the — > regeneration, the signal efficiencies for 
the default probability (about 0.1%) and about 0.2% or 0% 
(corresponding to ±100% error on the default probability) are 
compared and difference of the efficiencies is used as the sys- 
tematic error (less than 1% for all the searches). Total physics 
simulation errors are several percent in all four searches in 
SK-I, SK-II, and SK-III. In the background rate, major sys- 
tematic error sources used in the SK atmospheric neutrino data 
analyses fl7ll : neutrino flux, neutrino cross sections, and pion- 
nucleon interaction in oxygen, are taken into account and the 
total simulation errors are about 10-20%. 

As for the SK detector, uncertainty of the event selection 
cut value at each event selection step is considered as well as 
the energy scale uncertainty. Any detector uncertainty such as 
light attenuation and scattering in water is taken into account 
in the systematic error estimations which compare the cut ef- 
ficiencies between control sample data and MC. In the signal 
efficiency, the total detector errors are several percent in the 



searches. 
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FIG. 12: The K£ search; (upper plot) the signal efficiencies along event selections are shown in SK-I, SK-II, and SK-III from left to right, 
respectively ("FV" stands for the number of p — > p + K° MC events generated in true fiducial volume (22.5 kton). The events selected by the 
K® searches are excluded), (lower plot) the number of data candidates (black circle) and the expected atmospheric neutrino background rate 
(blue histogram) along event selections are shown in SK-I, SK-II, and SK-III from left to right, respectively ("FC" stands for the FCFV events 
described in Sec. IV). Only statistical errors are shown. There are no events remaining in the SK-II data after selection (D6). 



searches, the errors in SK-III are larger than those of SK-I 
by a factor of 2^3. The larger errors in SK-III are due to 
its poorer water quality which resulted in larger energy scale 
and ring counting systematics. The total detector error in the 
K*[ search (a few 10%) is mostly due to the vertex separation. 
The vertex separation error is estimated by changing the cut 
value corresponding to the difference between true and recon- 
struction described in Sec. V-A. In the background rates, the 
dominant error (about 20-60%) comes from statistics of the 
atmospheric neutrino events after all the event selections in 
the — >• 7t + 7r,~ Method 2 search. The other dominant error 
(about 20-40%) is from the vertex separation in the search. 



C. Comparison with Previous Result 

In this section we describe the differences between the 
present analysis and the previous result |5J]. For the estimation 
of the signal efficiency the previous analysis only considered 
the p — > /J + K® channel while the present work accounts for 
the entire p — > /j + K° channel including the K% contribution. 
This contribution coupled with improvements to the event re- 
construction has resulted in an efficiency increase from 5.4% 
to 7.0% in the K% -> K°K° search and from 7.0% to 10.6% in 
the K® — > n + K~ Method 1 search. Signal efficiencies and the 
background rates for the other searches are consistent within 
the systematic errors of the two analyses. 

Similar systematic error sources are incorporated in both 
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search method 


detector 


efficiency (%) 


background (%) 






total (physics, detector) 


total (physics, detector) 




SK-I 


7.5 (6.8, 3.2) 


14.2(11.0, 9.2) 




SK-II 


10.8 (6.8, 8.4) 


23.4(11.5, 20.4) 




SK-III 


10.0 (6.8, 7.4) 


19.1 (9.5, 16.5) 


K° -> n+n- Method 1 


SK-I 


9.8 (8.8, 4.2) 


15.8(15.2, 4.3) 




SK-II 


11.6 (8.8,7.6) 


22.8(16.5, 15.8) 




SK-III 


17.7 (8.9, 15.3) 


20.6 (16.2, 12.7) 


K% H> TC + JI~ Method 2 


SK-I 


7.5 (6.9, 2.8) 


63.7 (19.8, 60.5) 




SK-II 


7.1 (6.9, 1.8) 


34.4 (12.7, 32.0) 




SK-III 


7.9 (7.2, 3.4) 


24.6(11.6,21.7) 




SK-I 


17.2 (9.6, 14.3) 


31.9 (9.9, 30.3) 




SK-II 


19.2 (9.7, 16.6) 


38.9 (9.9, 37.6) 




SK-III 


18.4 (9.4, 15.8) 


48.4 (9.4, 47.5) 



TABLE II: Summary of systematic errors. The total systematic errors on the efficiency and background are shown in the third and fourth 
columns, respectively. The "physics" and "detector" in the parentheses represent the systematic errors from uncertainty of the physics simula- 
tions and the SK detector, respectively. 




FIG. 13: Typical event display of the real data candidate event in the 
search in SK-I. The cyan ring shows the reconstructed ring edge 
of the fj + candidate ring. The blue triangle and green circle show the 
reconstructed vertex position horizontally and vertically projected 
on the detector wall for the fj + and remaining particle candidates, 
respectively. The reconstructed vertex separation for this event is 
4.6 m. The reconstructed Cherenkov angle is 30.6° while the expec- 
tation from the jU-like momentum (291 MeV/c) is 37.6°. 



analyses. For instance, the error on the signal efficiencies 
is about 10% in both. However, the systematic error on the 
background rate in this analysis, ranging from 20% to 60% 
across the searches, is improved over the previous work's, 
which ranged from 40% to 80%. Most of this improvement 



is derived from the smaller statistical error on the atmospheric 
neutrino MC sample used in the present analysis. 

If the K9 search is excluded from the lifetime limit compu- 
tation and the data set restricted to the SK-I period, the result 
is slightly worse than the previous analysis, 1.1 xlO 33 years 
compared to 1.3 xlO 33 years. Despite the higher signal effi- 
ciencies and reduced systematic uncertainty on the expected 
background rates in this analysis, the increase in the number 
of observed candidate events in the — > 7T, + 7r.~ Method 1 
search degrades the lifetime limit. The three candidates found 
in the former analysis are also found here, but three additional 
events that were previously just outside of the signal region 
have now migrated in. Under the older reconstruction, their 
momenta fell slightly below the cut value of 250 MeV/c and 
thus failed to meet condition (B4) from Sec. VI-A. The im- 
proved reconstruction has increased their estimated momenta 
beyond this threshold. The p + candidate momentum recon- 
struction has been checked using the atmospheric neutrino 
MC and the shift in the reconstructed momentum relative to 
the true momentum, and the momentum resolution are both 
smaller in the present work. The new result represents an im- 
provement to the analysis. 



VII. CONCLUSION 

The SK p ->• p+K° decay search presented in has been 
updated to include improved event reconstruction algorithms 
and an additional 81.1 kton years of data. A search for pro- 
ton decay into has also been added. There is no sig- 
nificant excess of proton decay candidates found and the 
data remain consistent with the expected atmospheric neu- 
trino background rates for both the p — > p + K { ^ and p — »■ p + K® 
searches. Therefore a partial lifetime lower limit is set at 
x/B _^ + K o > 1.6 x 10 33 years at the 90% confidence level 
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using 172.8 ktonyears of data. The lifetime limits obtained 
by the K® and searches alone are 1.37 xlO 33 years and 
0.55 x 10 33 years, respectively. This result gives further con- 
straints on relevant SUSY GUT models. 
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